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Preface 


This  report  describes  procedures  followed  and  results  obtained  for  an 
evaluation  of  protection  provided  by  the  proposed  construction  of  a  deep-draft 
harbor  facility  in  the  ports  of  Ponce  and  Guayanilla,  Puerto  Rico.  Impacts  of  the 
proposed  construction  were  determined  by  conducting  numerical  simulations  of 
tidal  and  storm  surge  circulation  at  the  project  sites  using  with  and  without  the 
proposed  port  facilities.  Differences  in  computed  elevations  and  currents 
between  pre-  and  postconstruction  conditions  are  used  to  determine  project 
impacts.  Numerical  simulations  utilize  the  ADvanced  CIRCulation  (ADCIRC) 
long-wave  hydrodynamic  model.  This  study  was  performed  by  the  U.S.  Army 
Engineer  Research  and  Development  Center  (ERDC),  Coastal  and  Hydraulics 
Laboratory  (CHL),  at  the  Waterways  Experiment  Station  (WES)  in  Vicksburg, 
MS  for  the  U.S.  Army  Engineer  District,  Jacksonville  (CESAJ). 

The  investigation  reported  herein  was  conducted  by  Dr.  Norman  W. 
Scheffner,  research  hydraulic  engineer,  CHL.  Messrs  David  J.  Mark,  research 
hydraulic  engineer,  and  Fulton  C.  Carson,  computer  scientist,  CHL,  contributed 
significantly  to  the  study  by  developing  the  bathymetric  and  topographic 
database  and  installing  these  data  into  the  computational  grid.  The  study  was 
initiated  in  May  2001  and  completed  in  September  2001.  Mr.  David  Schmidt, 
CESAJ,  was  the  study  manager  and  point  of  contact  for  the  model  investigation. 
The  study  was  performed  under  the  general  supervision  of  Dr.  James  R.  Houston, 
former  Director,  and  Mr.  Thomas  W.  Richardon,  Director,  CHL. 

At  the  time  of  publication  of  this  report,  Dr.  James  R.  Houston  was  Director 
of  ERDC,  and  COL  John  W.  Morris  III,  EN,  was  Commander  and  Executive 
Director. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication, 
or  promotional  purposes.  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


1  Introduction 


The  U.S.  Army  Engineer  District,  Jacksonville,  requested  assistance  from  the 
U.S.  Army  Engineer  Research  and  Development  Center  (ERDC),  Coastal  and 
Hydraulics  Laboratory  (CHL),  in  evaluating  circulation  impacts  of  the  proposed 
development  of  a  megaport  facility  at  Ponce  and  Guayanilla  Harbors  along  the 
southern  coast  of  Puerto  Rico.  The  approach  taken  to  access  circulation  impacts 
utilizes  the  long-wave  hydrodynamic  model  Advanced  Circulation  Model  for 
Coastal  Ocean  Hydrodynamics  (ADCIRC)  to  compute  circulation  patterns  with 
and  without  the  proposed  expansion  projects.  Impacts  are  then  quantified  by 
computing  magnitude  changes  at  the  two  harbors.  Impacts  will  be  shown  for 
tidal  circulation  as  well  as  for  tropical  and  extratropical  storm  surge.  The 
following  sections  briefly  describe:  a)  the  ADCIRC  hydrodynamic  model; 
b)  verification  of  the  ADCIRC  model;  c)  storm  event  boundary  conditions; 
d)  project  impact;  and  e)  conclusions. 
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2  Hydrodynamic  Model 


In  the  following  two  subsections,  documentation  for  the  hydrodynamic 
model  is  first  given  to  provide  a  background  of  model  development  and 
demonstrate  the  applicability  of  the  ADCIRC  model  to  represent  the  study  area. 
This  section  is  followed  by  a  description  of  the  computational  grid  developed  for 
the  Puerto  Rico  study. 


Model  Documentation 

Water-surface  elevations  and  currents  for  both  tides  and  storm  events  are 
computed  through  applications  of  the  large-domain  long-wave  hydrodynamic 
ADCIRC  model  (Luettich,  Westerink,  and  Scheffner  1992).  The  ADCIRC 
model  is  an  unstructured  grid  finite-element  long-wave  model  developed  under 
the  U.S.  Army  Corps  of  Engineers  (USACE)  Dredging  Research  Program  (DRP) 
Griffis  et  al.  1995).  The  model  was  developed  as  a  family  of  two-  and  three- 
dimensional  codes  with  the  capability  of: 

a.  Simulating  tidal  circulation  and  storm  surge  propagation  over  large 
computational  domains  while  simultaneously  providing  high  resolution 
in  areas  of  complex  shoreline  and  bathymetry.  The  targeted  areas  of 
interest  include  continental  shelves,  nearshore  areas,  and  estuaries. 

b.  Representing  all  pertinent  physics  of  the  three-dimensional  equations  of 
motion.  These  include  tidal  potential,  Coriolis,  and  all  nonlinear  terms  of 
the  governing  equations. 

c.  Providing  accurate  and  efficient  computations  over  time  periods  ranging 
from  months  to  years. 

In  two  dimensions,  model  formulation  begins  with  the  depth-averaged 
shallow-water  equations  for  conservation  of  mass  and  momentum  subject  to 
incompressibility  and  hydrostatic  pressure  approximations.  The  Boussinesq 
approximation,  where  density  is  considered  constant  in  all  terms  but  the  gravity 
term  of  the  momentum  equation,  is  also  incorporated  into  the  model.  Using  the 
standard  quadratic  parameterization  for  bottom  stress  and  omitting  baroclinic 
terms  and  lateral  diffusion  and  dispersion,  the  following  set  of  conservation 
statements  in  primitive,  nonconservative  form  and  expressed  in  a  spherical 
coordinate  system  are  incorporated  in  the  model  (Flather  1988;  Kolar  et  al. 

1994): 
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where  t  represents  time,  X  and  <|>  are  degrees  longitude  (east  of  Greenwich  is 
taken  positive)  and  degrees  latitude  (north  of  the  equator  is  taken  positive),  x]  is 
the  free-surface  elevation  relative  to  the  geoid,  U  and  V  are  the  depth-averaged 
horizontal  velocities,  R  is  the  radius  of  the  Earth,  H=  C,  +  h  is  the  total  water 
column  depth,  h  is  the  bathymetric  depth  relative  to  the  geoid, /=  2Q  sin  9  is  the 
Coriolis  parameter,  Q  is  the  angular  speed  of  the  Earth,  ps  is  the  atmospheric 
pressure  at  the  free  surface,  g  is  the  acceleration  due  to  gravity,  r|  is  the  effective 
Newtonian  equilibrium  tide  potential,  po  is  the  reference  density  of  water,  ts>.  and 
are  the  applied  free-surface  stress,  and  t*  is  given  by  the  expression  C/JJ2  + 
V2)m  IH  where  C/ equals  the  bottom  friction  coefficient  which  can  be  specified  as 
either  linear  or  nonlinear  (Luettich,  Westerink,  and  Scheffiier  1992). 

The  momentum  equations  (Equations  2  and  3)  are  differentiated  with  respect 
to  X  and  t  and  substituted  into  the  time  differentiated  continuity  equation  (Equa¬ 
tion  1)  to  develop  the  following  Generalized  Wave  Continuity  Equation 
(GWCE): 
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The  ADCIRC-2DDI  model  solves  the  GWCE  (Equation  4)  in  conjunction 
with  the  primitive  momentum  equations  given  in  Equations  2  and  3. 

The  ADCIRC  model  solves  the  governing  equations  with  a  finite-element 
algorithm  over  arbitrary  bathymetry  encompassed  by  irregular  sea  and  shore 
boundaries.  This  algorithm  allows  for  flexible  spatial  discretizations  over  the 
entire  computational  domain  and  has  demonstrated  robust  stability  characteris¬ 
tics.  The  advantage  of  this  flexibility  in  developing  a  computational  grid  is  that 
larger  elements  can  be  specified  in  open-ocean  regions  where  less  resolution  is 
needed,  whereas  smaller  elements  can  be  applied  in  nearshore  and  estuarine  areas 
where  finer  resolution  is  required  to  resolve  hydrodynamic  details. 


Computational  Grid 

The  ADCIRC  model  has  been  applied  to  the  East  Coast,  Gulf  of  Mexico,  and 
Caribbean  Sea  domain  shown  in  Figure  1  to  develop  a  tidal  constituent  database 
(Westerink,  Luettich,  and  Scheffner  1993).  This  grid  was  used  as  a  basis  for  the 
present  Puerto  Rico  circulation  study.  A  subsection  of  the  grid  of  Figure  1  was 
first  developed  in  which  the  southern  boundary  was  maintained  at  the  coasts  of 
Venezuela  and  Columbia,  the  western  boundary  placed  midway  through  the 
island  of  Hispaniola,  and  the  northern  boundary  defined  as  an  arc  extending  north 
from  Hispaniola  and  eastward.  The  eastern  boundary  was  then  defined  by 
extending  the  arc  eastward  of  the  Lesser  Antilles  to  intersect  with  the  open  coast 
of  Surinam.  The  grid  and  boundaries  for  the  present  study  are  shown  in  Figure  2. 
Bathymetry  for  the  domain  was  initially  provided  from  the  National  Imagery  and 
Mapping  Agency  (NIMA  2000)  digital  database  of  bathymetry. 

Resolution  of  the  grid  in  the  study  areas  was  increased  to  provide  a  minimum 
grid  spacing  of  approximately  30  m.  This  resolution  allows  for  accurate  representa¬ 
tion  of  the  proposed  deep-draft  navigation  channel  and  docking  facility  at  Ponce  and 
Guayanilla.  The  increased  resolution  provided  for  the  harbors  is  shown  in  Figure  3. 
Existing  bathymetry  shown  for  Ponce  was  extracted  from  the  NIMA  digital  nautical 
chart  for  North  America  West.  Existing  bathymetry  for  Guayanilla  was  not  included 
in  the  NIMA  database  and  was  manually  entered  based  on  National  Oceanic  and 
Atmospheric  Administration’s  (NOAA)  National  Ocean  Survey  (NOS)  nautical 
chart  25681. 
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3  Verification  of  ADCIRC 
Model 


Verification  of  the  model  is  required  to  assure  that  grid  resolution,  bathymetry, 
and  boundary  conditions  are  adequately  prescribed  to  acceptably  reproduce  known 
or  observed  events.  Because  storm  surge  data  are  not  available  for  the  study  area, 
verification  is  limited  to  reproduction  of  tidal  surface  elevations  at  locations  within 
the  modeled  domain  shown  in  Figure  2  for  which  tidal  constituent  data  are  readily 
available. 

Tidal  propagation  is  simulated  within  ADCIRC  by  specifying  a  surface  eleva¬ 
tion  time  series  at  each  of  the  open-water  nodes  of  the  computational  grid  of  Fig¬ 
ure  2.  These  time  series  are  based  on  the  digital  tidal  constituent  database  of 
Le  Provost  et  al.  (1994)  for  the  M2,  S2,  N2,  Kj,  Oi,  Pi,  Qi,  and  K2  tidal  constitu¬ 
ents.  The  Le  Provost  database  has  been  successfully  used  in  all  recent  CHL  tidal 
studies  requiring  tidal  boundary  conditions  because  of  its  comprehensive  cover¬ 
age  of  the  world  and  its  consistent  accuracy  in  open  water.  The  ADCIRC  model 
was  run  with  these  boundary  conditions  for  a  15-day  simulation  period.  Com¬ 
puted  surface  elevation  time  series  were  archived  at  various  locations  in  the 
computational  domain  for  comparison  to  a  surface  elevation  time  series  recon¬ 
structed  from  existing  harmonic  analyses. 

Existing  harmonic  constituent  data  were  available  for  15  stations  within  the 
study  domain  from  the  NOAA  NOS  and  the  International  Hydrographic  Organi¬ 
zation  (IHO)  (1991).  Station  locations  are  listed  in  Table  1  and  shown  in 
Figure  4. 

Because  harmonic  data  are  not  available  at  Ponce  and  Guayanilla,  verification 
was  limited  to  demonstrating  that  the  model  acceptably  reproduces  tidal  elevations 
throughout  the  modeled  domain.  Therefore,  verification  is  limited  to  a  visual  com¬ 
parison  of  model  simulated  time  series  versus  reconstructed  time  series  based  on 
published  harmonic  constituents.  Figures  5  and  6  show  comparisons  of  simulated 
and  reconstructed  data  for  the  Puerto  Rico  coastal  stations  of  San  Juan  (sta  1)  and 
Magueyes  Island  (sta  12).  As  shown  in  the  figures,  the  simulation  provides 
sufficient  accuracy  to  satisfy  the  goals  of  the  present  study;  therefore,  it  is  concluded 
that  the  model  and  associated  grid  are  verified  for  tides.  Because  tides  and  surges 
are  both  examples  of  long  waves,  the  conclusion  that  the  model  is  verified  to  tides 
translates  into  the  assumption  that  the  model  can  acceptably  reproduce  tropical  and 
extratropical  storm  surge. 
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Table  1 

Station  Locations  \ 

Station 

East  Longitude 

North  Latatude 

Number 

Deg 

Deg 

Station  Name 

1 

-66.1167 

18.4617 

San  Juan,  La  Puntilla,  San  Juan  Bay,  PR 

2 

-64.7533 

17.6967 

Lime  Tree  Bay,  St  Croix,  VI 

3 

-64.8700 

18.3200 

Benner  Bay,  St.  Thomas,  VI 

4 

-64.9200 

18.3350 

Charlotte  Amalie,St  Thomas,  VI 

5 

-68.9333 

12.1000 

Willemstad,  Curacao 

6 

-66.9333 

10.6167 

La  Guaira,  Venezuela 

7 

-64.1667 

10.4500 

Cumana,  Venezuela 

8 

-61.5167 

10.6500 

Port  of  Spain,  Trinidad 

9 

-61.0000 

14.0167 

Castries,  St.  Lucia 

10 

-61.0500 

14.5833 

Fort-de-France,  Martinique 

11 

-64.9333 

18.3333 

St.  Thomas,  VI 

12 

-67.0500 

17.9667 

Magueyes  Island,  PR 

13 

-69.8833 

18.4667 

Ciudad,  Trujillo,  Dominican  Republic 

14 

-70.6833 

19.7500 

Puerto  Plato,  Dominican  Republic 

15 

-64.8833 

16.5333 

IAPSO  No.  30-1.3.2 
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Storm  Event  Boundary 
Conditions 


The  basis  of  determining  project  impacts  to  the  study  area  of  interest  was 
determined  to  be  magnitude  differences  in  surface  elevations  and  current 
magnitudes  as  a  result  of  the  proposed  project.  This  difference  was  to  be  defined 
for  both  tidal  circulation  and  tropical  and  extratropical  storm  surge.  Tidal 
boundary  conditions  are  based  on  the  tidal  constituent-based  boundary  conditions 
used  in  verification.  Storm  surges  were  selected  as  historic  events  that  impacted 
the  Island  of  Puerto  Rico.  The  following  two  subsections  describe  the  selection 
of  storm  events  to  be  used  in  the  storm  surge  comparison. 


Tropical  Storm  Surge 

The  hurricane  wind  field  model  used  in  conjunction  with  the  ADCIRC  model  is 
the  Hurricane  Planetary  Boundary  Layer  (PBL)  model  developed  by  Cardone 
(Cardone,  Greenwood,  and  Greenwood  1992).  This  model  simulates  hurricane¬ 
generated  wind  and  atmospheric  pressure  fields  by  solving  the  equations  of  hori¬ 
zontal  motion  which  have  been  vertically  averaged  through  the  depth  of  the  planet¬ 
ary  boundary  layer.  The  PBL  model  input  consists  of  histogram  and  snapshot  files 
which  define  the  hourly  location  in  latitude  and  longitude  of  the  eye  of  the  storm  and 
the  storm  intensity  parameters  specified  at  defined  times. 

Snapshot  and  histogram  files  are  computed  from  data  contained  in  the  NOAA 
National  Hurricane  Center's  HURricane  DATabase  (HURDAT)  of  tropical  storm 
events  (Jarvinen,  Neumann,  and  Davis  1988).  This  database  contains  descriptions  of 
all  hurricanes,  tropical  storms,  and  severe  tropical  depressions  which  have  impacted 
the  East  Coast,  Gulf  of  Mexico,  and  Caribbean  Sea  from  1886  to  present. 

A  tropical  storm  database  was  generated  during  the  DRP  (Scheffher  et  al.  1994) 
through  simulation  of  134  historically  based  storm  events  along  the  East  Coast,  Gulf 
of  Mexico,  and  Caribbean  Sea  based  on  the  HURDAT  database.  For  486  discrete 
locations  along  the  U.S.  coast,  peak  storm  surge  values  corresponding  to  storm 
events  which  produced  a  surge  of  at  least  0.305  m  were  archived  and  indexed 
according  to  event,  location,  and  surge  magnitude.  This  indexed  database,  as  well  as 
the  most  current  version  of  the  HURDAT  data,  were  used  to  select  Hurricane 
Georges,  which  crossed  Puerto  Rico  on  22  September  1998,  as  the  event  for 
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determining  hurricane  impact  to  the  study  area.  A  plot  of  the  track  of  Hurricane 
Georges  is  shown  in  Figure  7. 


Extratropical  Storm  Surge 

In  an  approach  similar  to  that  of  the  tropical  event  database  previously 
described,  an  extratropical  storm  event  database  was  generated  within  the  DRP. 

This  database  was  constructed  by  driving  the  ADCIRC  model  with  wind  fields 
extracted  from  the  U.S.  Navy  Fleet  Numerical  Meteorology  and  Oceanography 
Center's  database  of  winds  for  the  1 6-year  winter  storm  period  (defined  as 
September  through  March)  of  1977  through  1993  (77-78,  78-79,  etc.).  These  wind- 
field  data  are  provided  at  a  6-hr  temporal  interval  on  a  2.5°  latitude  and  longitude 
spatial  grid.  The  extratropical  storm  database  consists  of  a  7-month  surface 
elevation  and  current  hydrograph  at  each  of  the  486  stations  previously  described. 

Time  series  plots  corresponding  to  an  archived  station  near  the  center  of  the 
study  area  were  analyzed  for  peak  values  in  surface  elevation  and  current  magnitude. 
Each  time  series  represented  surge  with  no  tide.  The  time  series  plot  corresponding 
to  DRP  sta  682,  located  offshore  of  Ponce  and  Guaynilla,  was  generated  for  each  of 
the  1 6  storm  year  periods.  One  period  of  extratropical  surge  was  found  to  occur 
during  the  extratropical  storm  year  1992-1993  shown  in  Figure  8.  In  this  figure,  day 
1  refers  to  1  September  1982  such  that  the  storm  event  of  day  195  refers  to  15  March 
1993. 
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5  Project  Impact 


Project  impacts  were  defined  to  be  changes  in  surface  elevation  and/or  depth- 
averaged  currents  as  a  result  of  the  construction  of  the  proposed  harbor  facilities.  In 
the  following  two  sections,  the  specifics  of  both  proposed  projects  are  summarized. 
This  summary  is  followed  by  results  of  the  numerical  simulations  for  with  and 
without  project  construction.  For  each  project  scenario,  contours  of  maximum 
surface  elevation  and  maximum  current  magnitude  and  vector  differences  are  used  to 
quantify  project  impacts. 


Ponce  Harbor 

The  proposed  deep-draft  harbor  complex  at  Ponce  Harbor  involves  the 
deepening  of  the  navigation  channel  to  a  minimum  of  13.7  m  and  development  of  a 
new  dock  facility  in  the  harbor.  This  new  dockage  will  require  the  reclaiming  of 
approximately  60  acres  of  water  to  the  immediate  north  of  Punta  Penoncillo.  The 
area  of  reclaimed  water  is  shown  in  the  project  footprint  superimposed  on  the 
topographical  map  in  Figure  9.  The  model  representation  of  the  reclaimed  feature  is 
shown  in  Figure  10.  It  is  assumed  that  the  reclaimed  area  is  not  subject  to  wetting 
and  drying  but  will  remain  a  permanent  land  feature.  Minimum  depths  in  the 
navigation  channel  are  set  to  15.0  m. 


Guayanilla  Harbor 

The  proposed  deep-draft  harbor  complex  in  Guayanilla  Harbor  involves 
construction  of  an  extensive  harbor  facility  on  landfill  just  north  of  the  end  of  Punta 
Guayanilla  as  shown  in  Figure  11.  The  off-loading  facility  adjoining  the  landfill  will 
be  on  piers  and  not  represent  an  obstruction  of  flow.  Model  representation  of  the 
reclaimed  feature  is  shown  in  Figure  12.  It  is  assumed  that  the  reclaimed  area  is  a 
permanent  land  feature  and  not  subject  to  wetting  and  drying.  The  off-loading  pier 
is  located  in  deep  water  so  that  no  dredging  is  required. 


Simulation  Results  -  Tides 

Harbor  construction  impact  on  tidal  circulation  was  based  on  an  analysis  of  a 
28-day  tidal  month  simulation  over  the  full  domain.  The  tide  was  based  on  the  eight 
constitutuents  used  in  model  verification  and  did  not  include  winds.  Resulting  time 
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series  of  tidal  elevations  and  currents  were  archived  at  the  15  reference  locations 
shown  in  Figures  13  for  Ponce  and  Guayanilla. 

Tides  along  the  southern  coast  of  Puerto  Rico  are  primarily  diurnal  with  one 
high  tide  per  day  as  shown  in  the  time  series  for  Magueyes  Island  in  Figure  6. 
Because  all  reference  stations  are  in  unprotected  open  water  and  the  two  harbors  are 
adjacent  to  each  other,  reference  station  tides  are  essentially  identical.  For  example, 
the  28-day  surface  elevation  time  series  for  stas  8  and  13  are  shown  in  Figure  14. 
Note  in  the  figure  that  the  first  2  days  show  a  reduced  amplitude  due  to  starting  the 
model  from  mean  sea  level  (MSL)  and  ramping  the  tide  to  full  amplitude  in  3  days. 
Therefore  tides  representative  of  the  study  area  should  be  taken  from  days  3-28.  The 
tidal  current  magnitude  is  shown  in  Figure  15  for  the  corresponding  28  days. 

Tides  in  both  Ponce  and  Guayanilla  Harbors  are  small  with  a  maximum  spring 
tide  amplitude  on  the  order  of  0.15  m  and  currents  less  than  1  cm/sec.  However,  in 
order  to  quantify  construction  impacts  over  the  entire  harbor  area,  time  series  of 
surface  elevation  and  currents  are  archived  at  all  computational  nodes  of  the  model 
to  facilitate  generation  of  maximum  magnitude  difference  contour  maps  and 
difference  vector  maps.  Although  tides  are  not  affected  much  by  local  changes  to  the 
shoreline  configuration  or  bottom  bathymetry,  the  impact  of  the  two  harbor  projects 
are  shown  in  Figure  15  for  tidal  surface  elevation  extreme  (maximum  or  minimum) 
differences  and  Figure  16  for  maximum  tidal  current  magnitude  vector  differences. 
The  time  period  for  computed  maximum  differences  was  selected  as  the  spring  tide 
period  of  days  12.5  through  17.5  of  Figure  14.  In  Figure  15,  differences  represent 
the  maximum  elevation  before  construction  less  the  maximum  elevation  after 
construction  of  the  landfills  (i.e.,  the  maximum  difference  occurring  during  the  5-day 
spring  tide  of  Figure  14).  Similarly,  the  vectors  and  vector  magnitudes  shown  in 
Figure  1 6  represent  the  preconstruction  current  vector  less  the  postconstruction 
vector  at  the  time  of  greatest  difference.  Therefore,  the  magnitude  of  the  vector 
represents  the  reduction  or  increase  in  currents  as  a  result  of  construction  of  the 
proposed  harbor  facilities. 

Inspection  of  Figures  16-18  shows  that  the  proposed  harbors  generate  minimal 
impact  on  tidal  circulation.  For  example,  maximum  elevation  changes  at  either 
Ponce  or  Guayanilla  are  less  than  0.01  cm  with  corresponding  current  magnitude 
differences  of  less  than  1 .0  cm/sec.  These  small  impacts  are  to  be  expected, 
however,  due  to  the  slow  change  of  tide  for  the  diurnal  signal  and  the  fact  that  wind 
forcing  is  not  included  in  the  tidal  simulations.  Wind  impacts  are  addressed  in  the 
following  section  in  the  tropical  storm  event  simulation  of  Hurricane  Georges. 


Simulation  Results  -  Tropical 

Impacts  of  the  proposed  projects  can  best  be  evaluated  by  looking  at  extreme 
events.  For  this  reason.  Hurricane  Georges  was  selected  as  the  test  tropical  event. 
Hurricane  Georges  passed  over  the  Island  of  Puerto  Rico  on  22  September  1998 
generating  a  positive  storm  surge  in  both  Ponce  and  Guayanilla  Harbors  due  to  the 
wind-forcing  of  water  into  the  harbors  from  the  southwest  quadrant  of  the  storm  as  it 
passed  north  of  the  study  areas.  Storm  surge  elevations  at  open  water  stas  8  and  13 
are  shown  in  Figure  19. 
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Storm  surges  in  the  Caribbean  Sea  are  generally  small  as  shown  in  the  surge 
values  of  Figure  19.  This  is  due  to  the  fact  that  surges  in  the  Caribbean  propagate 
from  deep  water  immediately  into  shallow  water  with  no  broad  continental  shelf  to 
force  shoaling.  Therefore,  large  surges  characteristic  of  the  East  Coast  and  Gulf 
Coast  of  the  United  States  do  not  occur  in  the  Caribbean  Sea.  However,  wind-driven 
circulation  can  be  significant  in  shallow  harbors  and  protected  areas,  therefore, 
circulation  impacts  of  the  proposed  projects  are  used  to  evaluate  project  impact. 

Figures  20  and  21  show  wind-driven  circulation  vectors  in  Ponce  Harbor  with 
and  without  the  proposed  harbor  expansion  computed  at  the  time  of  the  peak  of  the 
storm  (day  3.28  of  Figure  19).  Similarly,  Figures  22  and  23  show  circulation  vectors 
in  Guayanilla  Harbor  with  and  without  the  proposed  expansion.  Figures  20-24  also 
indicate  surface  elevation  contours  that  show  maximum  surge  elevations  of  less  than 
0.4  m.  Figures  24  and  25  show  maximum  current  difference  vectors  computed 
during  the  life  of  the  event  superimposed  on  maximum  current  magnitude  difference 
contours  for  Ponce  and  Guayanilla  Harbors.  Figures  23  and  24  demonstrate 
maximum  impacts  of  the  proposed  harbors  and  indicate  current  reductions  to  be  on 
the  order  of  up  to  1 .0  m/sec.  Maximum  reductions  are  shown  to  be  at  Guayanilla 
Harbor  where  the  preconstruction  surge  circulation  vectors  of  Figure  21  are  shown 
to  be  deflected  offshore  in  Figure  22.  However,  impacts  of  the  proposed  harbor 
expansion  are  limited  to  the  close  proximity  of  the  project. 


Simulation  Results  -  Extratropical 

The  extratropical  event  of  March  1993  previously  described  was  selected  to 
represent  a  winter  storm  event  impact  to  Ponce  and  Guayanilla  Harbors.  Because 
extratropical  event  winds  are  much  less  intense  than  tropical  events,  the  resulting 
surge  is  found  to  be  considerably  smaller  than  for  a  tropical  event.  For  example. 
Figure  26  shows  the  surface  elevations  for  preconstruction  existing  conditions  at 
stas  8  and  13  for  the  10  days  of  simulation  (10-20  March  1993).  As  shown,  the 
surge  elevations  are  less  than  10  cm.  The  corresponding  wind-driven  currents  near 
the  proposed  harbors  at  Ponce  and  Guayanilla  were  each  on  the  order  of  0.2  m/sec 
as  shown  in  Figures  27  and  28  for  day  8.125  of  Figure  26. 

Computed  maximum  differences  in  current  magnitudes  and  vectors  for  Ponce 
and  Guayanilla  are  shown  in  Figures  29  and  30  to  be  less  than  0.20  m/sec.  Currents 
and  current  vector  differences  between  pre-and  postconstruction  for  the  extratropical 
event  are  much  less  than  the  corresponding  computations  for  the  tropical  event; 
therefore,  it  is  concluded  that  maximum  impacts  resulting  from  construction  of  the 
proposed  harbor  complex  should  be  based  on  the  tropical  event  simulations. 
Conclusions  based  on  this  analysis  are  presented  in  the  following  section. 
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6  Conclusions 


ERDC,  CHL  was  requested  by  the  Jacksonville  District  to  determine  how 
proposed  deep-draft  facility  expansions  of  the  harbors  of  Ponce  and  Guayanilla 
along  the  south  coast  of  Puerto  Rico  may  impact  the  hydrodynamics  of  the 
surrounding  coastal  environment  The  ADCIRC  long-wave  hydrodynamic 
model  was  used  to  estimate  tidal  propagation  and  storm  surge  in  each  of  the 
harbors  for  both  preconstruction  existing  conditions  and  postconstruction  future 
condition  deep-draft  harbors.  Impacts  of  the  proposed  harbor  expansions  were 
determined  by  computing  surface  elevation  and  current  differences  between  pre- 
and  postconstruction.  Events  selected  for  simulation  were  a  typical  28-day  lunar 
month,  Hurricane  Georges,  and  an  extratropical  event  of  March  1993. 

Both  Ponce  and  Guayanilla  Harbors  have  southern  exposure  that  is  partially 
protected  by  offshore  islands  and  shallow  regions.  Offshore  bathymetry  drops  to 
over  600  m  within  5-10  km  of  the  entrance  to  either  harbor.  Because  of  the 
exposure  and  lack  of  an  offshore  shelf,  tides  and  storm  surges  do  not  become 
well  developed,  but  remain  small.  For  example,  spring  tides  (without  wind)  in 
either  harbor  are  less  than  0.2  m  in  amplitude  and  the  maximum  storm  surge  for 
Hurricane  Georges  was  less  than  0.4  m.  Therefore,  surface  elevation  impacts  of 
the  proposed  harbors  are  small.  For  example,  the  maximum  project  impact  on 
surface  elevation  for  Hurricane  Georges  at  Guayanilla  Harbor  was  less  than 
0.05  m. 

Wind-driven  currents  within  the  bay  represents  the  most  potential  long-wave 
threat  to  the  coastal  infrastructure  resulting  from  the  passage  of  a  tropical  or 
extratropical  event.  For  this  reason,  change  in  hurricane  surge  currents  as  a  result 
of  construction  of  the  harbor  expansions  was  identified  as  the  best  measure  of 
construction  impact.  Pre-and  postconstruction  differences  in  current  magnitudes 
were  computed  as  a  means  of  demonstrating  reductions  or  increases  in  current  as 
a  result  of  the  proposed  landfill. 

Computations  indicate  that  preconstruction  depth-averaged  currents  in  the 
vicinity  of  the  Ponce  landfill  were  on  the  order  of  0.5  m/sec  during  peak 
positive  surge.  For  Guayanilla  Harbor,  peak  currents  at  peak  surge  were  also  on 
the  order  of  0.75  m/sec.  Computed  maximum  differences  in  current  at  any  time 
during  the  storm,  i.e.,  peak  positive  surge  to  peak  drawdown,  were  0.6  -  0.7 
m/sec  at  Ponce  and  0.9  -  1.0  m/sec  at  Guayanilla.  The  values  represent  decreases 
in  velocity  as  a  result  of  the  harbor  expansions  because  the  constructed  landfill 
deflects  the  current  offshore  of  the  landfill.  These  reductions  were  shown  to  be 
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localized  to  the  point  that  other  regions  within  the  harbor  were  unaffected  by  the 
construction. 

Conclusions  of  this  study  are  twofold.  First,  simulations  show  that  the 
harbors  of  Ponce  and  Guayanilla  do  not  experience  large  tides  or  tropical  and 
extratropical  storm  surges.  Secondly,  pre-  and  postproject  simulations  of  a 
severe  tropical  event  show  that  projects  impacts  are  small  (less  than  1.0  m/sec) 
localized  reductions  in  storm  circulation  currents  which  do  not  affect  regions 
further  than  approximately  1  km  from  the  proposed  projects. 
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Figure  2.  Computational  domain  of  Puerto  Rico  study 


a.  Ponce 


b.  Guayanilla 

Figure  3.  Grid  resolution  and  existing  bathymetry 
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Figure  4.  Station  locations  within  model  domain 


Figure  5.  Comparison  of  ADCIRC  and  IHO/NOS  reconstructed  tide  at  San  Juan 
(sta  1 ) 
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Figure  6.  Comparison  of  ADCIRC  and  IHO/NOS  reconstructed  tide  at  Magueyes 
Island  (sta12) 


Figure  7.  T rack  of  Hurricane  Georges 


Figure  8.  Surge  values  for  extratropical  storm  year  1 992-1 993  at  DRP  sta  682 


Figure  10.  Model  representation  of  proposed  Ponce  Harbor 


Figure  1 1 .  Topographic  map  of  Guayanilla  with  project  footprint 


Figure  12.  Model  representation  of  proposed  Guayanilla  Harbor 


Figure  1 3.  Reference  stations  for  Ponce  and  Guayanilla  Harbors 
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Figure  1 5.  Tidal  current  magnitude  preconstruction  time  series  for  28  days, 
stas  8  and  13 


Figure  16.  Maximum  tidal  elevations  differences  during  days  12.5-17.5 
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Figure  1 8.  Maximum  tidal  current  differences  during  days  1 2.5-1 7.5  for  Guayanilla 
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Figure  1 9.  Storm  surge  elevation  for  Hurricane  Georges  at  stas  8  and  1 3 


Figure  20.  Pre-expansion  circulation  vectors  at  Ponce  for  Hurricane  Georges  at  day  3.28 


Figure  21 .  Post-expansion  circulation  vectors  at  Ponce  for  Hurricane  Georges  at  day  3.28 


Figure  22.  Pre-expansion  circulation  vectors  at  Guayanilla  for  Hurricane  Georges  at  day  3.28 


Figure  23.  Post-expansion  circulation  vectors  at  Guayanilla  for  Hurricane  Georges  at  day  3.28 


Figure  24.  Maximum  current  difference  at  Ponce  for  Hurricane  Georges 


Figure  25.  Maximum  current  difference  at  Guayanilla  for  Hurricane  Georges 
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Figure  27.  Existing  condition  extratropical  storm  circulation  vectors  at  Ponce  at  day  8.125 


Figure  29.  Maximum  extratropical  current  difference  at  Ponce 


Figure  30.  Maximum  extratropical  current  difference  at  Guayanilla 
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